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DISPERSIVE TREATMENT OF /sT ^ tttt 
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We discuss a new method to treat the K -kit amplitude dispersively, taking into full account 
the effects of final state interatcions. Our approach is based on a set of dispersion relations for 
the K — > TTTT amplitude, in which the weak Hamiltonian carries momentum. In these dispersion 
relations two subtraction constants have to be introduced, whereby one can be related via a soft 
pion theorem to the K ^ n amplitude. The second is presently unknown, and we use lowest 
order Chiral Perturbation Theory for a first guess of it's value. We emphasize the advantage of 
combining this approach with lattice input which could provide the two subtraction constants 
with sufficient accuracy. 



1 Introduction 



Lattice QCD would in principle be the appropriate tool for the calculation of the K vrvr ampli- 
tude, accounting fully for its nonperturbative nature. In practice, an explicit calculation of 
K TTTT is not yet feasible, although there was some progress in this direction rece ntlylll. The 
major obstacle for direct lattice calculations were already identified some time ago by Maiani 
and Testa in a no-go theorem, stating that for decays into two or more particles, their interaction 
with each other makes a direct calculation on the lattice impossible@. The standard workaround 
for K — > TTTT is the calculation of the unphysical K ^ n matrix element, and the use of a relation 
valid at lowest order in chiral perturbation theory (CHPT), to relate it to the K — > vrvr matrix 
element!. Since lowest order SU{3) CHPT is known to be accurate only at the 30% level, the 
step from — > vrto K — > vrvr induces large uncertainties. Unfortunately, the usage of a one loop 
CHPT relation is not possible since this would imply the use of several low energy constants 
which are not available!. 

A crucial point for the calculation of the K vrvr matrix element is the inclusion of final 
state interaction (FSI), whose importance in the context of e^e was pointed out by Pallante 
and PichBQ, following the ideas outlined in a paper of TruongB, who showed that the inclusion 



of FSI for K —5- vrvr decays yields an enhancement of the / = amphtude, pointing in the right 
direction concerning the AI = 1/2 rule . These FSI are totally neglected if one relates the 
K ^ TT amplitude with the K — > vrvr amplitude with the help of the tree level CHPT relation. In 
order to to write down a dispersion relation for the K — > vrvr amplitude, Pallante and Pichi use 
an offshell kaon field, which can be defined in infinitely many ways, introducing also ambiguities 
in the final numerical result for the K vrvr amplitude. A detailed discussion of this point may 
be found in ref. @. 



2 Dispersive treatment with momentum carrying Hamiltonian 

One can avoid the problems related to the use of an offshell kaon field by allowing the weak 
Hamiltonian to carry momentum; a procedure which has been suggested in ref. u. We will sketch 
how the method works: Define the amplitude 

l=0{APiMP2)\n]l\0)\K{qi))=:T+{s,t,u) , (1) 

with the Mandelstam variables s = {pi+P2)'^, t = (qi—pi)'^, u = ^2)^1 related by s + 
t + u = 2M^ + + §2 J where q2 is the momentum carried by the weak Hamiltonian. The 
physical decay amplitude is obtained by setting q2 = {s = M'^, t = u = M^). To describe 
a function of three variables dispersively would be very complicated. The problem can be 
simplified considerably if we neglect the contribution of the imaginary parts of D waves and 
higher. In this approximation, the amplitude decomposes into several functions depending each 
only on one of the three Mandelstam variables: 
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where A = — M^. Mq[s) corresponds to / = 5- wave in the s channel, whereas in the t 
channel Nq and A'^i denote the 1 = 1/2 5' and P wave and Rq the 1 = 3/2 5' wave. 

The dispersion relation of the full amplitude is converted into a set of coupled dispersion 
relations of functions of a single variable, which can be solved numerically. For instance, for 
Mo(s), giving the major contribution in the final result, we define the right-hand cut: 

discMo(s) = sin 5^{s)e-< \Mo{s) + Mo{s) 

and get the dispersive representation: 



Mo(s) = O[;(s,so) i« + ^(s-so 
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with two subtraction constants (SC) a and b. 



The Omnes function 0,q{s, sq) is defined to be: 



i2o(s, soj = exp < / as 

vr J4M2 
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Mo(s) is an angular average of Mo(s). Likewise, one can define the same quantities for the 
remaining functions A'o,! and iioi but no more new SC's have to be implemented!. 

The number of SC's which have to be introduced is crucial. In order to solve the dispersion 
relation and calculate the K tttt amplitude, we have to provide two SC's as input. One of 
those can be linked to the K ^ n amplitude: a soft pion theorem relates the amplitude at the 
soft pion point {s = u = M^,t = Mj^) to the K —>■ t: amplitude up to order corrections: 

-^%^ = a.^^OiMi), (3) 

where N = Nq{M]^) + 2Rq{M]^). Notice that although the process involves a Kaon, the sym- 
metry argument leading to the above relation is based on SU{2), and suffers therefore only from 
0{M^) corrections. With the help of this relation, the first subtraction constant can in principle 
be provided by a lattice calculation. The problem is the second subtraction constant b; it is 
related to the derivative in s of the amplitude r+ at the soft pion point: 

b = (s, e-s,m2) 1,^^,2 + ... 

In ref. §1 a Ward identity which relates this derivative to a -fC ^ tt matrix element is derived. 
The calculation of this matrix element would provide b. 

Another option to get b is to calculate the amplitude (|^) for the following unphysical kine- 
matical values: 

s = AMl t = u=^-Ml , (4) 

where the two pions are produced at rest0. This specialkinematic configuration does not get 
into conflict with the no-go theorem of Maiani and TestaH. 

In the absence of a value for b, one can illustrate the numerical results by fixing it at a value 
and then varying it within a fairly wide range. For its central value one can use lowest order 
CHPT: 



3a 



Ml - M2 



(l + X + 0{M%)) . (5) 



The size of the correction X is at the moment unknown, but nothing protects it of beeing of the 
order of M|. 

The numerical analysis is shown in Fig.|l], where T'^{s, Mj^ — M) as a function of the incoming 
Kaon momentum squared, s, is plotted. For the uncertainty X coming from next to leading 
order CHPT , a rather wide range X = ±30% has been chosen. Comparison with the lowest 
order CHPT formula, also plotted in Fig.|l], shows that large corrections have to be expected 
due to the Omnes factor, if we neglect next to leading order CHPT effects (X=0). If we vary X 
in the above given range, we see that we enhance the effect for positive X and decrease it for 
negative X. In the case X = —0.3 the Omnes-corrected curve differs only little from the lowest 
order CHPT one. 
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Figure 1: The function \T^ {s,t,u)\ plotted vs. E = y/s along the hue of constant u = M^: the result of our 
numerical study for different values of X are compared to tree level CHPT. 

This analysis shows again that in order to get a reliable number for the K vrvr matrix 
element in this framework, it is crucial to obtain exact values for the two subtraction constants. 
The accuracy of the final result is essentially limited by these, since the other potential source 
for uncertainties, the phase shifts needed as input for the Omnes functions, are known to a 
rather high precision. 
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